source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
energy = wdi[is.na(wdi$energy_pc_gj) == F]
energy_interp = energy[, list(interpolate(year, energy_pc_gj)) , by = code3]
observation = lengths(energy_interp$V1)
g = do.call(rbind, energy_interp$V1)
View(g)
test = rep(energy_interp$code3, observation)
View(energy_interp)
observation = lengths(energy_interp$V1)
energy_interp$code3
test = rep(energy_interp$code3, observation)
test = rep(energy_interp$code3, each = observation)
sum(observation)
sum(observation/2)
observation = nrows(energy_interp$V1)
observation = lengths(energy_interp$V1) / 2
test = rep(energy_interp$code3, each = observation)
test = rep(energy_interp$code3, observation)
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
test = data.table(codes, data)
View(test)
plot(test$year, test$value)
View(wdi)
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
energy = wdi[! is.na(energy_pc_gj)]
View(energy)
energy = wdi[! is.na(energy_pc_gj)][,1:4]
View(energy)
names(oecd)
oecd_names = names(oecd)
oecd_names = names(oecd)
oecd_names[1]
oecd_names[1] == "code3"
oecd_names[1] = "code3"
oecd_names
oecd_names[12] = "year"
oecd_names
names(oecd) = oecd_names
oecd_energy = merge(oecd, energy, by = c("code3", "year"))
View(oecd_energy)
plot(oecd_energy$energy_pc_gj, oecd_energy$Value)
View(oecd_energy)
un_wiid_names = names(un_wiid)
oecd_names[2] = "code3"
names(un_wiid_names) = un_wiid_names
un_wiid_eergy = merge(un_wiid, energy, by = c("code3", "year"))
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
un_wiid_energy = merge(un_wiid, energy, by = c("code3", "year"))
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
plot(un_wiid_energy$energy_pc_gj, un_wiid$Gini)
plot(un_wiid_energy$energy_pc_gj, un_wiid_energy$Gini)
plot(un_wiid_energy$energy_pc_gj, un_wiid_energy$Gini, log ="xy")
plot(un_wiid_energy$energy_pc_gj, un_wiid_energy$Gini, log ="x")
View(un_wiid_energy)
boxplot(un_wiid_energy$Gini)
boxplot(un_wiid_energy$Gini[un_wiid_energy$energy_pc_gj >= 50])
year
test = un_wiid_energy[, mean(Gini), by = .(code3, year)]
test = un_wiid_energy[, mean(Gini), by = .(energy_pc_gj)]
plot(test)
plot(test, log = "x")
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
hist(un_wiid$Gini)
merge = rbind(
data.table(country = oecd_energy$Country,
code = oecd_energy$code3,
year = oecd_energy$year,
energy = oecd_energy$energy_pc_gj,
gini = oecd_energy$Value
)
)
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
View(merge)
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
View(merge)
View(un_wiid)
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
plot(merge$energy, merge$gini)
plot(merge$energy, merge$gini, log = "x")
source('~/Desktop/origin_inequality/Supplementary Material/Empirical Data/Energy Gini Merge/energy_gini.R')
